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Abstract 



Antineutrinos born in the U and Th decay chains inside the Earth 
■ ( "Geoneutrinos" ) carry out information on the amount and distribu- 

(—] ' tion of radiogenic heat sources, which is of fundamental importance for 

Qh! geophysics. Models of the Earth distribute U and Th masses mainly 

Q^' between the continental crust and the lower mantle. It has been much 

^ • discussed recently that a number of detectors stationed at appropriate 

geographical sites can separate the crust and mantle contributions. In 
present work we analyze directional separation of iJ^ signals arriving 
^ ■ from the crust and the lower mantle with only one detector. We find 

that with a ~30-kton liquid scintillation antineutrino spectrometer 
using i7e + P — > n + e"*" detection reaction and positron and neu- 
tron coordinates reconstruction techniques the U and Th distribution 
model can roughly be tested. We also consider detector calibration us- 
ing ~1 MCi commercially available ^''Sr-^'^Y beta source which emits 
i7(, {Eyrnax = 2.28 MeV) in the geoneutrino energy range. 

Introduction 

The Earth emanates about 40 TW of heat flow coming from the interior [1]. 
It is beheved that considerable part of this heat originates from radioactive 
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decay of Uranium, Thorium and Potassium hidden in the Earth. Antineutri- 

nos z7g born in the radioactive decays ( "Geoneutrinos" ) bring information on 
the abundances and radiogenic heat sources inside the Earth, which are of a 
key importance for understanding of the formation and subsequent evolution 
of our planet. 

The Geoneutrino concept is now 45 years old. A short overview of this 

concept can be found in our paper [2]. Only recently it was recognized that 
Geoneutrinos emitted in the Uranium and Thorium decay chains can be 
detected in large volume liquid scintillation underground spectrometers [3,4] 
using the inverse beta decay process 

+ P ^ + (1) 

as detection reaction. Actually KamLAND GoUaboration in Japan has al- 
ready demonstrated revolutionary progress in detection technique and has 
reported that a few geoneutrino events have been observed [5]. BOREXINO 
GoUaboration at LNGS, Italy, is planning experiments in the same field. 

U and Th abundances and their distributions in the Earth are not known, 
except for a thin layer near the surface where direct sampling is possible. 
Some information is obtained from analysis of volcanic outflows. Orthodox 
model of the Earth (sometimes labeled as BSE model) and geochemichal 
arguments distribute U and Th masses mainly between the continental crust 
and the lower mantle and thus the geoneutrino fiux is expected to vary with 
geographical position of the observational point, decreasing from continental 
to oceanic sites (Table 1). 



Table 1: Expected rates of U + Th geoneutrino events at various sites (ex- 
position 10^^ H • year*, no oscillations) 



Site 


Rothschild et al.** 


Raghavan et al. 


Fiorentini et al. 




[3] 


[4] 


[6,7] 


Himalaya 


65 




112 


Baksan 






91 


Gran Sasso 


53 


la: 134; lb: 50 


71 


Kamioka 


48 


la: 75; lb: 50 


61 


Hawaii 


27 




22 



* 1160 ton of CH2 contain 10 hydrogen atoms. 
** Calculated using Uf. fluxes from Ref. [3]. 
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As can be seen in Table 1, at continental sites dominate geoneutrinos 
born in the crust, at oceanic sites antineutrino flux is a few times lower. 

Geoneutrino data from several experiments around the world can sepa- 
rate crust and mantle contributions [3,4,6,7], and thus to confirm or discard 
predictions of the Earth's model. Experiment of this type has been proposed 
for Baksan Neutrino Observatory in North Caucasus [2]. 

The CHOOZ experiment in 1999 y demonstrated that the use of reaction 
(1) to detect low energy {E < 9 MeV) reactor z7g in liquid scintillation detec- 
tor provides determination of z7g incoming direction in case of (practically) 
one point like i7e source [8]. 

In this paper we study separation capability of geoneutrino {E < 3.27 
MeV) signals coming from different directions (from the crust and the mantle) 
using one ~30-kton target mass liquid scintillation detector. At Baksan 
sample of ~4000 geoneutrino events can be accumulated in ~5 year (live 
time) data taking. We conclude with detector calibration procedure using 
90gj._90Y \)Qta, source, which emits z7g {E^j^ax — 2.28 MeV) in the geoneutrino 
energy range. 



1 Neutrino Directions 

We consider detection in (CH2)n liquid scintillator with a density of 
0.8 g/cm^. Geoneutrinos detection signature is delayed coincidence between 
the prompt positron signal and signal from the neutron capture gammas. 
Positron energy spectrum boosted by absorption of two 0.511 MeV anni- 
hilation quanta is shown in Fig. 1. For each event positron and neutron 
capture coordinates (Rgj and R„i) are reconstructed using scintillation light 
collection pattern and (or) distribution of light arriving times detected in 
photomultipliers . 

Information on neutrino incoming direction(s) is obtained from studies of 
the reconstructed positron- neutron position vectors Reni- 

Reni Rei Rni (2) 

Below we present a number of Monte Carlo experiments, which include 
modeling of Geoneutrinos coming from mantle, from continental crust and 
from both of these sources at a time. To check our Monte Carlo simulations 
we start with geoneutrino parallel beam coming from a remote point-like 
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source already considered at CHOOZ for the reactor u^s [8]. For reaction (1) 
differential cross section and kinematics we use results presented in [9,10]. 

1.1 Parallel geoneutrino beam 

Low energy neutrons emitted in reaction (1) (their energies are in the 1—15 
keV range) are strongly concentrated around incoming geoneutrino direction, 
which is assumed as Z axis (Fig. 2). In subsequent collisions the memory 
of this direction is partially conserved and neutrons are displaced from the 
reaction (1) point in the positive Z direction. Formation of this displacement 
in the first few collisions can be traced step by step in Fig. 3a,b. After first 
~ 8 collisions the memory is completely lost and neutrons are slowing down 
and diffuse symmetrically around the displaced center. Slowed down neutrons 
are captured mostly in hydrogen with mean capture time r ~ 200 /is. The 
energy independent neutron - ^^C scattering cross section is 4.8 b, neutron - 
proton scattering cross section is 20 b for neutron energies grater than 1 eV, 
for lower energies due to hydrogen binding in the molecule the cross section 
increases gradually and reaches 65 b at thermal energy of 1/40 eV. Mean 
square radius < >^/^ of the neutron capture distances from the (displaced) 
center is calculated to be 6.1 cm. The average neutron displacement dz in Z 
direction is found to equal 1.72 cm: 

dz = 1.72 cm, (r^)^/^ = g.l cm (3) 

The main factor, which limits the sensitivity of the method, is the accu- 
racy of neutron - positron position reconstruction procedure. We assume that 
vector (2) projections are Gaussian-distributed around their mean values: 

(Ren,) = ± ajN'/\ {Keny) = ± a./iV^/^, 

(Ren.) = 1.72 cm ± ajN^/^, a^ = ay = a, = a (4) 

where N is the number of detected reaction (1) events and cr is the Gaussian 
dispersion. 

Assuming now cr = 20 cm and N — 2500 neutrino events (as in the 
CHOOZ experiment) we find (Renz) = 1.7 ± 0.4 cm, and thus neutron dis- 
placement can be detected at four standard deviation level. In polar coordi- 
nates the most probable direction of the average vector (Rem) is ^ = 0; the 
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uncertainty of this value is S9 ~ 2^^^a/1.7 • N^^^ ~ 18° in agreement with 
CHOOZ [8]. 

Another way of presenting resuhs is to build, in polar coordinates, the an- 
gular distribution W{cos9) of reconstructed vectors Rem around i7e incoming 
direction, which can be presented in the form: 

W^(cos 9) ^ Const ■{l + kz-cos9) (5) 

In the case considered here kz ~ {dz/cr). (In Eqs. (5) a higher order 
term a • cos'^9 with a ~ {dz/crY is neglected). We generate 10^ z7g events 
and find true angular distribution: W(cos9) — Const ■ {1 -\- kz ■ cos 9) , with 
kz = 0.138 (solid line in Fig. 4). Next we generate 2500 events and find 
the uncertainty d^^ = ±0.03, represented by shaded area in Fig. 4. Thus 
the angular distribution (5) differs from isotropy at four standard deviation 
level. We conclude that both ways of analysis (displacement and angular 
distribution) are equivalent. 

1.2 Geoneutrinos from the Lower Mantle 

Under the crust is the Earth's mantle, which is 2900 km thick. Its lower part 
(the lower mantle) is believed to contain as much U and Th as the Earth's 
crust [6,7]. Geoneutrino flux has finite vertical component and is assumed 
to be symmetric azimuthally (Fig. 5). Again we generate 10^ Monte Carlo 
events and find true value of neutron cloud displacement diM = 1.20 cm in 
the vertical direction. Thus 

(Renx) = (Reny) = ± a/N^^^, 

dLM = (Ren.) = 1-20 Cm ± (t/ATV^. (g) 

With 4000 detected neutrino events c^lm = 1-2 ± 0.32 cm can be found. 

1.3 Geoneutrinos from the Crust 

Here we consider hypothetical case where continental crust source forms a 
uniform 6000 km diameter and 40 km thick circular region centered around 
the geoneutrino detector. In this model average horizontal components of 
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reconstructed positron-neutron vector Re„ vanish, the vertical component is 
small: 



The displacement for crust geoneutrinos dcr is 4 times smaller than that 
for the neutrinos from the lower mantle and cannot be obsereved with the 
number of detected geoneutrinos N = 4000 considered here. 

1.4 Geoneutrinos from the Crust and the Lower Mantle 

We consider detector installed at a continental site (e.g. at Baksan Observa- 
tory). The average displacement of neutron cloud in the vertical direction is 
given by the expression: 

C^LM+Cr = (R-enz) = aLM^^LM + (1 " ^Lm) " C?Cr ± Ct/N'^/^, (8) 

where o^lm ~ 1-2 cm, dcr ~ 0.29 cm and olm = Fi^m/{Fi^m + -^Cr) is the 
lower mantle fraction in the total geoneutrino incoming flux; N — 4000, the 
maximal achievable events sample considered here, a jN^I'^ = 0.32 cm. 

Dependence (8) is shown in Fig. 6. One can sec that separation method 
considered here is not very sensitive. Only sufficiently large displacements 
c^LM+Cr, larger than 1 cm if found, can indicate contradiction to the pre- 
dictions of the orthodox Earth's model. In case experiment favors lower dis- 
placements, and thus indicates low contribution of the mantle geoneutrino 
flux, the dominant role of the crust geoneutrinos predicted by the model can 
roughly be confirmed. Only with much larger number of collected events 
(A'' ~ 2 • 10'^) and therefore with much larger detector more definite conclu- 
sions could be reached. 

1.5 Detector Calibrations 

Detection of small displacements discussed above requires adequate cahbra- 
tion procedures. While usual method, based on inserting neutron and gamma 
sources into the fiducial volume can and should be exploited, use of suf- 
ficiently strong z7g source is highly desirable. We propose to consider for 



(Ren.) = (Reny) = ±(7/7V^/^ 




(7) 
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calibration purposes movable ~1 MCi ^°Sr-^°Y beta source. ^°Sr (T1/2 = 
28.6 yr) decays to the ground state of ^^Y{E^raax = 2.28 MeV, T1/2 = 64 h), 
which with 99.99% probability populates the ground state of the stable ^^Zr 
nucleus. If installed at the distance of 30 m from the 30 kton detector center, 
the source can generate about 2 • 10^ events of reaction (1) per year [11]. 

Two circumstances make this source attractive. First, it will irradiate 
the detector with z7g flux of known intensity, known energy spectrum in the 
geoneutrino energy range and of known angular structure and, second, the 
sources are produced commercially and used to supply heat for Radioisotope 
Thermoelectric Generators (RTGs) (see Fig. 7). 

We note that proposed cahbration method could also be used in other 
low energy i7g experiments employing large liquid scintillation detectors. 

2 Discussions and Conclusion 

In this paper we present a first attempt to analyze directions of a multidirec- 
tional low energy z7g flux. In this attempt only two geoneutrino sources have 
been taken into account: the continental crust and lower mantle. We haven't 
analyzed perturbations due to possible fluctuations of U and Th concentra- 
tions in the detector's immediate vicinity. Clearly more work is needed to 
come to more accurate results. 

At this preliminary stage of analysis we can summarize the results as 
following: 

• Present understanding of the total amount of radiogenic sources and 
their distribution in the Earth's reservoirs is based on a shaky ground 
of cosmogenical and geochemical arguments with an obvious deficit of 
direct experimental evidence [1, 6, 7, see also discussion in Ref. [2]]. 
Information obtained with one 30-kton target mass detector using di- 
rectional separation of incoming geoneutrino fiux is useful but limited: 
it can give only some indications against the orthodox Earth's model 
predictions, or provide it's rough confirmation. More definite informa- 
tion could be obtained only with much larger detector. To avoid mis- 
understanding we note that the Georeactor hypothesis can be tested 
at BNO with ~l-kton antineutrino detector [13]. 

• With 30-kton scintillation spectrometer additional information on geoneu- 
trinos can be obtained: The total geoneutrino fiux can be measured 
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within ~2% accuracy; accurate spectroscopy of positron spectrum (see 
Fig. 1) will provide important information on global U to Th concen- 
tration ratio [12]. "Geoneutrino" is a part of future studies at BNO of 
low energy z7g fluxes of natural origin. Among other physical targets 
we mention: Galactic Supernova explosion and estimation of frequency 
of gravitational collapses in the Universe by detection of isotropic flux 
of Relict antineutrino (for more points see Ref. [2] and project LENA 
[14]). 
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Figure 1. Expected geoneutrino induced positron spectrum at Baksan [2], 

32 

(exposition 10 H year, no oscillations). Solid line - U + Th, 
1 - Th, dashed line is estimated signal from power reactors. 
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Figure 2. Angular distribution of reaction (1) neutrons relative to the incident 
geoneutrino direction weighted with the reaction cross section 
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Figure 3. Monte Carlo stimulation of the positron- neutron displacement 

growth in a few first neutron collisions with scintillator atoms, 
(a) shows position of the first collision, <Z> indicates the average 
coordinate of the 1st- collision, (b) The growth of <Z> with 
the number of collisions. Geoneutrino beam is along Z axis. 
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Figure 4. Angular distribution of the projections of positron-neutron reconstructed 
vector Ren along the incident direction. Solid line is "true" distribution 

found with 10^ Monte Carlo events, shaded area represents uncertainty of 
the slope for 2500 events sample. 




Figure 5. Internal structure of Earth (schematic). Lines show trajectories 
from the lower mantle to the detector (D). 
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Figure 6. Vertical displacement of neutron cloud d|_|^^Q|. vs relative contribution 
of lower mantle to the total (crust + lower mantle) geoneutrino 
flux (solid line). Shaded is the (68% CL) uncertainty region, dark grey area 
between vertical lines represents model' prediction for detector installed 
in BNO. 



